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ABSTRACT 

Wc obtained optical long-slit spectra of four planetary nebulae (PNe) with low- 
ionization pair of knots, namely He 1-1, IC 2149, KjPn 8 and NGC 7662. 

These data allow us to derive the physical parameters and excitation of the pairs 
of knots, and those of higher ionization inner components of the nebulae, separately. 

Our results are as follows. 1) The electron temperatures of the knots are within 
the range 9 500 to 14 500 K, similar to the temperatures of the higher ionization 
rims/shells. 2) Typical knots' densities are 500 to 2 000 cm~'^. 3) Empirical densities 
of the inner rims/shells are higher than those of the pairs of knots, by up to a factor of 
10. Theoretical predictions, at variance with the empirical results, suggest that knots 
should be denser than the inner regions, by at least a factor of 10. 4) Empirical and 
theoretical density contrasts can be reconciled if we assume that at least 90% of the 
knots ' gas is neutral (likely composed of dust and molecules). 5) By using [Raga et al.l 
(|2008l ) shock modeling and diagnostic diagrams appropriated for spatially resolved 
PNe, we suggest that high-velocity shocked knots traveling in the photoionized outer 
regions of PNe can explain the emission of the pairs of knots analysed in this paper. 

Key v^fords: (ISM): Planetary nebulae: individual: He 1-1, IC 2149, KjPn 8, 
NGC 7662. Planetary nebulae: structure. 



1 INTRODUCTION 

The small-scale low-ionization structures (LIS) of planetary 
nebulae were classified in terms of their morphology and 
kinematics by Gongalves, Corradi & Mampaso (2001) as 
knots, jets and jetlikeS systems, eithe r in pairs or isol ated. 
From this, and subsequent works (e.g. lGoncalve3l2004l . and 
references therein) it has became apparent that: i) around 
10% of the Galactic PNe are known to possess LIS; ii) they 
are indistinctly spread among all the PNe morphological 
classes; iii) 50% of these PNe have highly coUimated, high- 
velocity jets, and/or high- velocity pairs of knots (FL IERs; 
fast, low-ionization emission regions. iBalick et al.lll993 ): and 
iv) most of them are mainly photoionized. 

However, a number of key questions remain open: are 
there significant density contrasts between LIS and the 
higher ionization main components (rims, attached shells 



^ These are features whose morphology resemble highly super- 
sonic jets, but that share the expansion of the nebular compo- 
nents in which they are embedded. 



and detached haloes) of the nebula? What are the typical 
electron densities of the different types of LIS? 

We are involved in a long-term project for characteriz- 
ing the small-scale LIS of PNe since they can tell us much 
about the formation and evolution of PNe. The main issues 
we can access by studying the LIS are: 

(i) the collimation processes of LIS in regard to those pro- 
cesses responsible for the shape of the PN itself, or the mass- 
loss processes during the AGB and post-AGB phases; 

(ii) the effect of the ionization front on the fossil AGB 
features; and 

(iii) the role of disks and magnetic fields in shaping the 
highly coUimated outfiows (jets) in PNe. 

Our motivation for the present analysis is to go a step 
further on the study of the different LIS analyzing their den- 
sities. LIS can move through their environments either with 
high or low velocity. In the first case (the so-called FLIERS), 
radial velocities of 24 — 200 km s~^ with re spect to the 
PN m ain components are typically measured (|Balick et al.l 
I1993D . whereas in the second instance (the SLOWERs, slow 
moving low-ionization emission regions; see |Perinottdl2000l ) 
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Table 1. Log of the INT+IDS observations 



PN Name and P.A. 


Obs. Date 


Exposure (s) 






3x 


He 1-1 - 315° 


Aug 31 


300, 2400 


IC 2149 - 70° 


Sep 01 


300, 900, 1800 


KjPn 8 - 98° 


Aug 30 


1200 


KjPn 8 - 120° 


Aug 30 


300, 1200 


NGC 7662 - 175° 


Sep 05 


60, 300 


NGC 7662 - 248° 


Sep 04 


30, 60, 300 



the knots do not show pecuhar velocities and share the typi- 
cal 30 kms~^ (for elliptical shells) to 100 kms~^ (for bipolar 
lobes) expansion velocities of their host PNe. A straightfor- 
ward question arises from this: how is the density contrast 
between the LIS and rims, shells or haloes correlated with 
the different velocity regimes? 

In this paper we analyse 4 PNe that contain pairs of 
knots, either of high- or low- velocities. 

We present our spectroscopic data and treatment, in 
Section 2. The line fluxes obtained for as many as possible 
large- and small-scale components of the four PNe, and the 
derived physical properties for each of the components are 
shown in Sections 3 to 6. Finally, Section 7, is dedicated 
to the discussion of the results in terms of the match with 
the theoretical model predictions, and to our concluding re- 
marks. 



2 OBSERVATION AND DATA TREATMENT 

A 120 sec exposure Hq + [N II] (6568/95A) image of 
He 1-1 (Figure [1} was retrieved from the IPHAS database 
(INT/WFC Photometric Ha Survey of the Northern Galac- 
tic Plane: |http://www.i phas. org ). We also retrieved from the 
Hubble Space Telescope (HST) Archive two WFPC2 images 
of IC 2149: the V-band (filter 555W, 2 sec exp. time) and 
the [O III] (F502N, 100 sec) images, obtained on 1995 (pro- 
posal 6119), see Figure The [N ll] (F658N, 400 sec) image 
of NGC 7662 (Perinotto et al. 2004; Figure [5]), was also re- 
trieved from the HST archive, and was obtained on 1996. 
The image of KjPn 8 was obtained with the 2.5-m Nordic 
Optical Telescope (NOT) at the Observatorio del Roque de 
los Muchachos (European Northern Observatory, La Palma, 
Spain), with ALFOSC camera, in June 26, 2002. The image 
we show in Figure |4] is the result of one exposure of 1 800 s 
through the narrow- band filter [N ll]6584A. 

Spectra of the 4 PNe were obtained on August 31 and 
September 1, 4 and 5 of 2001, at the 2.5 m Isaac New- 
ton Telescope (INT) at the Observatorio del Roque de los 
Muchachos using the Intermediate Dispersion Spectrograph 
(IDS). The 235 mm camera and the R300V grating were 
used, providing a spectral coverage from 3650 to 7000 A with 
a spectral reciprocal dispersion of 3.3 A pixel" The spatial 
scale of the instrument was 0'.'70 pixel" ^, with the TEK5 
CCD. Seeing varied from 0'.'9 and I'.'l. The slit width and 
length were l'.'5 and 4 ', respectively. These data were taken 
with the slits positioned through the centre of the nebulae at 
the position angles (P.A.), and with exposure times given in 
Table 1. Each of the exposures was taken 3 times, as listed 
in the table. The longer exposures were used to measure the 



fluxes of most emission lines. However, for almost all the 
PNe, the [O iii]AA 4959,5007 A, the [N ii]AA 6548,6583 A 
doublet as well as the Hq and even the H/3 lines were satu- 
rated at the brightest features in the longer exposures, thus, 
the shorter exposures where used in these cases. Differently 
of the way we took long-slit spectra of the other PNe, in the 
case of KjPn 8 the slit was centred on the position of the 
central star, and then it was shifted in order to cover each of 
the knots. Compared to the weather conditions during the 
observation of the other three PNe, that of the KjPn 8 was 
poor (clouds appeared during the integration). 

During the night, bias frames, twilight and tungsten 
flat-field exposures, wavelength calibrations, and exposures 
of standard stars (BD +332642, Cyg OB2 No. 9, HD 19445, 
and BD -1-254655) were obtained. Spectra were reduced 
following the IRAF instructions for long-slit spectra, be- 
ing bias-subtracted, flat-fielded, combined in order to im- 
prove the signal-to-noise ratio (S/N) and eliminate cos- 
mic rays, wavelength-calibrated, and sky-subtracted. Fi- 
nally, they were flux-calibrated using the above mentioned 
standard stars and the mean atmospheric extinction curve 
for La Palma. 



3 HE 1-1 

PN G055. 34-02.7 (He 1-1) is a point-symmetric PN with a 
pair of low-ionization knots. See in Figure [1] the extensions 
of its bright Rim (or a combination of a rim and a shell), 
the large-scale component that dominates the emission of 
the nebula, and the SE and NW strings of knots, which 
correspond to the small-scale faint structures seen in the [N 
ll] NOT image. 

In the uppermost part of Table 2 we list selected line 
fluxes of the Rim, the knots and the entire portion of the 
nebula covered by the slit. These include all the electron 
density and temperature diagnostic emission lines covered 
by the IDS spectra. Absolute H/3 fluxes, Fh/3, integrated 
along the slit for each nebular component are also given 
in Table 2. The Balmer lines ratios (Hq/H/3, H7/H/3 and 
H5/H/3) were then used to derive the Ch/3, the logarithmic 
ratio between observed and dereddened H/3 fluxes. For the 
derivation of ch/3, we assumed Te = 10* K and the densities 
given by the [S 11] ratio of each region . Theoretical Balmer 
line ratios fr om [Osterbrock fc Ferlandl (l2006t) a nd the red- 
dening law of lCardeUi. Clavton fc MathisI l| 19891 ) were used. 
The weighted average of Ch/3, per PN component, is the one 
given in the table, ch/3 is constant along the slit, and the 
ch/3 = 1.49±0.16 for the e ntire nebula is l o wer th an the value 
previously published by iTvlenda et al.l (|l992l ). The latter 
authors found a range of values for ch/3, which goes from 
1.9 to 2.2. Note, however, that their values are based only 
on the Ha/H/3 ratio. Fluxes were then dereddened using the 
derived ch/3. The dereddened fluxes, or intensities, are not 
shown in the table. Instead, flgures in the table correspond 
to the observed fluxes. 

The middle part of Table 2 shows the errors on the 
fluxes. They were calculated taking into account the statis- 
tical errors in the measurements, as well as systematic errors 
of the flux calibrations, background determination, and sky 
subtraction. 

In the lowermost part of Table 2 electron densities and 
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Table 2. He 1-1: Observed line fluxes, flux errors, cjj/s, Ne and Te 



Figure 1. He 1-1: IPHAS H«-(-[N ll] image. The 3 structures un- 
der analysis are indicated in the image. The size of the extraction 
windows in the spectra includes 4'.' 2 each side from the centre for 
the Rim, and extends from -|- 4'.'9 to 7", and from -4'.'9 to -7'.'7 
from the centre, for the SE and the NW knots respectively. The 
extraction window for the entire portion of the nebula covered by 
the slit, NEB, extends along 14^7. 



Line ID 


Rim 


Knots 


NEB 








NWk 




[Sll] 4072.0 


6.232 


- 


- 


6.464 


H6 4101.8 


15.57 






15.97 


H7 4340.5 


29.69 


- 


- 


29.13 


[Olll] 4363.2 


11.36 


16.89 


14.61 


13.04 


H/3 4861.3 


100.0 


100.0 


100.0 


100.0 


[Olll] 4958.9 


543.4 


587.7 


601.7 


566.8 


[Olii] 5006.86 


1715. 


1773. 


1808. 


1774. 


[CLill] 5517.7 


3.321 


22.56 


15.57 


4.152 


[CLin] 5537.9 


3.041 


18.23 


13.04 


3.640 


[Nil] 5754.6 


11.90 


23.46 


24.87 


13.83 


[Nil] 6548.0 


319.2 


480.6 


413.9 


335.1 


Ha 6562.8 


993.9 


977.6 


981.2 


1021. 


[Nil] 6583.4 


1034. 


1436. 


1321. 


1075. 


[Sll] 6716.5 


129.3 


114.4 


124.4 


132.3 


[Sll] 6730.8 


169.6 


114.7 


139.2 


172.1 


Ph/3 " 


21.62 


0.4565 


0.6247 


22.37 


CH/3 


1.46±0.20 


1.44±0.32 


1.45±0.33 


1.49±0.16 






Percentage error 


s in line fluxes 




(0.01-0.05)H/3 


35 


52 


60 


31 


(0.05-0.15)H/3 


15 


25 


27 


13 


(0.15-0.30)H/3 


11 


20 


20 


11 


(0.30-2.0)H/3 


10 


17 


16.6 


8.5 


(2.0-5.0)H/3 


9.5 


15 


16 


6.5 


(5.0-10.0)H/3 


9 


14.5 


16 


6.5 


> lOH/3 


9 


14.5 


16 


6.6 




Electron 


Densities (em~'' 


) and Temperatures (K) 


Ne[S II] 


1600±230 


600±140 


900±210 


1560±190 


Nejci III] 


1700±830 


800±230 


1000±330 


1360±580 


Te [O III] 


12500±2060 


14400±3320 


13450±4110 


13060±1580 


Te [N 11] 


10800±1800 


12900±2980 


14000±3320 


11400±1590 


Te[S II] 


8750±1450 






9300±1350 



In units of 10 



within the errors, whereas Ne of the Rim is higher than that 
of the SE and NW knots, by factors of ~2.7 and ~1.8, re- 
spectively. There is a good agreement between the two tem- 
perature diagnostics obtained for the three PN components. 

Figure [2] shows that Te at the centre could be some- 
what lower than at the knots, but the variation is within 
the errors. 

Summarising, the pair of knots in He 1-1 are substan- 
tially less dense than the Rim, whereas no strong evidence 
exist for significant temperature variations from the Rim to 
the knots. 



temperatures are presented. At least two different measures 
of Te and Ne are available for each PN component. They 
represent regions of low (Ne[S ll], Te[N ll] and Te[S ll]) and 
medium excitation (Ne[Cl ill] and Te[0 ill]). 

The densi t y of the nebula can be compared with 
ISamland etall (l992h . They find Ne[S Il] = 1.19xl0^ cm~^ 
in moderate agreement with ours. Figures for Te and Ne 
of the knots we re previously publish ed only in a proceeding 
contribution bv lBem'tez et al.l (|2002l '): 11 300 K for the tem- 
perature of the NW knot, and densities of 300 cm~'^ and 
600 cm~'^, for the SE and NW knots respectively. These val- 
ues hardly agree with ours. Note that they did not derive, 
but assumed the 10* K of the SE knot. 

Figure [2] (top left panel) and Table 2 shows that both 
density measurements in each region of the nebula are equal 



4 IC 2149 

PN G166. 1-1-10. 4 has an apparent shape that does not easily 
compare with the classification bins (round, elliptical, bipo- 
lar or quadrupolar, irregular and point-symmetric) usually 
adopted for P Ne. Kinematic mode l ling suggests that it is 
a bipolar PN ([Vazquez et al.l |2002| : iFeibelman et all Il994l : 
IZhang fc Kwoklll998l '). As shown in Figure^ IC 2149 is com- 
posed by a bright higher ionization emission zone, the Core, 
and by a pair of lower excitation knots, the North-Eastern 
one being much brighter than its counterpart to the South- 
West. In fact only the two former stru ctures can be identifie d 
from Figure [1 Please refer to Fig.l of lVazguez et all l|2002l ). 
in which not only a higher excitation image is show, but also 
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NGC 7662; P.A.= 175 
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Figure 2. Electron densities and temperatures as a function of the distance to the centre of the nebula. The different structures under 
analysis are shown as symbols that correspond to the centre of each structure, as defined in the images (Figures 1, 3, 4, and 5). Filled 
(Ne[S ll]) and open (Ne[Cl III]) triangles represent densities, while filled (Te[N ll]) and open (Te[0 III]) circles, as well as the filled box 
(Te[S ll]) represent electron temperatures. Symbols are plotted slightly displaced in distance in order to avoid overlapping. Horizontal 
lines represent densities and temperatures averaged along the slit (NEB): continuous line for Ne[S ll], large-small dashed line for Ne[Cl 
III], small dashed line for Te[0 III], medium dashed line for Te[N ll], and large dashed line for Te[S ll]. TOP LEFT: He 1-1 - from left to 
right the NW knot, the Rim and the SE Imot, respectively, as in Figure [T] BOTTOM LEFT: IC 2149 - From left to right, the NE knot, 
the Core and the SW knot as in Figure |3] TOP RIGHT: NGC 7662 - P.A.=175°; five zones of the Inner rim and Outer shell, as shown 
in Figure|5]and Table 5. BOTTOM RIGHT:NGC 7662 - P.A.=248°; further five zones of the Inner rim and FLIERs, following Figure|5] 
and Table 5. 



the lower excitation [N ll] and [S ll] are presented. In their 
[N ll] image the SW LIS is also clearly seen. 

Table 3 contains the spectroscopic results for IC 2149, 
obtained from the INT-I-IDS spectra at P.A.=70° (observed 
and absolute H/3 fluxes, flux errors, ch/3, Ne and Te) for the 
Core and pair of LIS of the nebula. 

The Ha/H/3 ratio was used to obtain the ch/3 of the 
NEk, SWk and NEB regions. Owing to t he presence of 
Balmer line absorption in the central zones (jVazguez et al.l 



l2002l l. the Ch/3 of the Core was adopted as the average 
from the knots' values. The value measured for the nebula 
ch3=0-12±0.01 is in agreement with other p ublished values 
(iFeibelman et al ] | 19941 : ICiarduUo et"al]|l999h but it is much 
smalle r than the value (0.41) obtained by IVazguez et al.l 
ll2002t) for the Core and the knots. 

The gas physical conditions we derived to this PN (Ta- 
ble 3 and left bottom p anel of Fig. 2) compar e relatively well 
with those derived bv lVazguez et aEI l|2002f ) for the North- 
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P.A.=70 




Table 3. IC 2149: Observed line fluxes, flux errors, ch/j, Ne and Te 



Figure 3. IC 2149: HST 555W (V) image. Since this filter is 
centred in [O III], only the Core is clearly identified in this image. 
The low-ionization NEk and its fainter SWk counterpart are hard 
to see in t his image, but are eas ily identified in the [N ll] image 
m Fig.l of lVazQ uez et al. The spectroscopic sizes of the 

Core, NE and SW knots are, respectively: 2''l each side from the 
centre; - 2'.'8 to - 8'.'4; and + 2'.'8 to 6'.'3. The entire PN, NEB, is 
14'.'7 in size. 



Eastern knot and for the NEB. The exception is the Ne[S 
ll], for which they obtained values of 10 000 and 9 600 cm""^ 
while we found 5 350±1 550 and 4 650±660 cm-^ for the 
NEk and NEB, respectively. We cannot understand the rea- 
son of this discrepancies, since they do not depend on the 
very different ch/3 applied to correct the fluxes. On the 
other hand, the Te that would be affected by our different 
choices in terms of cha, are in reasonable agreement with 
theirs f see I Vazquez et al.|[200^ Table 1). If only NEB is con- 
cerned, our resu lts also agree with the results obtained by 
iFeibelman et al.l (|l994)- 

As [CI III] AA5517,5537A were not detected, only Ne[S 
ll] could be derived for the Core. Figure [5] shows that the 
two density diagnostics for each of the knots are such that 
SWk (Ne[S ll]=2 050±1 130) is significantly less dense than 
NEk, whose density is similar to that in the Core (Ne[S 
ii]=5 350±1 550 and 6 050±1 450, respectively). 

In spite of an apparent increase of Te[0 ill] and Te[N 
ll] from the North-Eastern to the South- Western side of the 
nebula, variations are within the errors. 



5 KjPn 8 

PN G112. 5-00.1 (KjPn 8) is a huge PN (size ~ 14 'x 
4 '; Lopez, Vazquez & Rodriguez 1995) with multiple po- 
lar ejections. A number of papers deal with the character- 
istics of this PN, including: i) narrow- band imaging of it s 
many structures in the difi'erent scales, iLopez et al.l (|l995h: 
ii) pro per motions of its latest bipolar ejections, MeaburnI 
(| 19971 ): iii) kinematics of the large bipolar lobes and of the 



Line ID 


Core 
ore 


Knots 


NEB 






NEk 


SWk 




[Sll] 4072.0 


- 


5.743 


5.625 


3.011 


H6 4101.8 


15.75 


26.99 


36.27 


24.14 


H7 4340.5 


40.46 


43.41 


61.99 


46.12 


[Olll] 4363.2 


3.838 


2.457 


5.333 


3.850 


H/3 4861.3 


100.0 


100.0 


100.0 


100.0 


[Olll] 4958.9 


180.1 


124.5 


140.7 


150.9 


[Olii] 5006.86 


540.2 


378.4 


408.4 


454.0 


[OljIIIJ OOi i . I 




0.367 


0.545 


0.386 


[cijinj 000 /.y 




0.516 


0.481 


0.440 


[Nil] 5754.6 


1. 140 


1.633 


1.227 


1.299 


[Nil] 6548.0 


15.03 


30.27 


18.55 


22.06 


na 0062.0 


320.6 


321.3 


277.8 


315.1 


[Nil] 6583.4 


43.07 


91.54 


51.81 


65.23 


[Sll] 6716.5 


1.288 


3.201 


2.418 


2.261 


[Sll] 6730.8 


2.312 


5.607 


3.385 


3.860 


Ph/3 " 


11.17 


12.78 


4.549 


31.76 


CH/3 


0.07'' 


0.14±0.015 


0.00 


0.12±0.008 






Percentage error 


s in line fluxes 




{0.01-0.05)H/3 


17 


24 


39 


10 


{0.05-0.15)H/3 


10 


17 


24 


5.5 


{0.15-0.30)H/3 


7 


11 


16 


4.5 


{0.30-2.0)H/3 


6.5 


8 


10 


4.5 


{2.0-5.0)H/3 


5.5 


7 


8.5 


4.5 


{5.0-10.0)H/3 


4.5 


6.5 


7 


4 


> 10 H/3 


4.5 


6 


6.5 


4 




Electron 


Densities (cm^^ 


) and Temperatures (K) 


Ne[S II] 


6050±1450 


5350±1550 


2050±1130 


4650±660 


Ne[CI III] 




7600±2750 


1500±870 


4000±680 


Te[0 III] 


10350±1760 


10150±2560 


12700±3150 


11050±1170 


Te [N 11] 


12300±2150 


10450±2600 


12350±4850 


11000±1170 



In units of 10 



values is the average between the knots' values; see text. 



pair of knots. iLopez et al] |l993); and iv) temperature, den- 
sity and chemical abundances of the inner nebula, Vazquez, 
Kingsburgh & Lopez (1998). At much smaller scales, a ring 
that seems to coUimate the many structures (with its CO 
and H2 counterpart) has been detected giving support to 
the idea that KjPn 8 had two PNe-like events in its history 
(|L6pez et al.l (I2OO0I ')'). However, Ne and Te of its two pairs 
of low-ionization knots were not properly measured so far, 
due to their extreme faintness. 

In our spectra only one of the knots of each pair could be 
measured with a good S/N. At variance with the other PNe 
in our sample, no analysis of the internal or entire nebula 
will be given, because of the huge extension of the PN. Our 
results for the two SE knots of KjPn 8 are given in Table 4. 

When comparing our ch;3 values with those in the lit- 
erature, particular care should be taken because they usu- 
ally refer to different nebu lar regions, and redd ening may 
vary along this large nebula. [Vazguez et al.l l|l998l ) measured 
very different ch/3 values for the three knots which they la- 
belled Al, A3 and Bl; ch/j = 1-49, 0.96 and 0.56, respec- 
tively. Our Ch/3 for SEki2o (their Al) is therefore 3 times 
smaller than theirs, whereas for SEk98 (their knot Bl) both 
determinations agree very well, ch,3=0.56 and 0.52±0.10, 
respectively. The reason for the discrepant ch,3 in SEki2o 
is not clear at present. As for the density, they reported 
Ne[S ll]=100 cm"^ for both SEkgs and SEki2o, while we de- 
rived substantially larger values Ne[S ll]=600±140 ([CI lll] = 
450±270) and 600±90 for SEk98 and SEki2o, respectively. 
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Figure 4. KjPn 8: NOT [N ll] image. The two P.A. of the struc- 
tures under analysis are indicated in the image. The spectroscopic 
sizes of the two knots, SEkgg is 14" and that of SEki2o is lO'.'S. 



Table 4. KjPnS: Observed line fluxes, flux errors, ch/3, Ne and Te 



Line ID 


South Eastern Knots 




SEki20 


SEkgs 


[Sll] 4072.0 


10.78 


8.869 


4101.8 


23.35 


20.69 


H7 4340.5 


41.41 


42.54 


[Olll] 4363.2 


1.149 


2.466 


H/3 4861.3 


100.0 


100.0 


[Olll] 4958.9 


100.5 


101.8 


[Olll] 5006.86 


307.9 


307.4 


[CLlll] 5517.7 




3.000 


[CLlll] 5537.9 




2.280 


[Nil] 5754.6 


9.351 


12.95 


[Nil] 6548.0 


392.8 


398.0 


Ha 6562.8 


438.8 


453.3 


[Nil] 6583.4 


1191. 


1207. 


[Sll] 6716.5 


99.79 


98.28 


[Sll] 6730.8 


100.4 


98.19 


Fh/3 " 


2.468 


1.221 




o.48±a.a7 


0.52±0.1O 




Pcrccnta^ 


■•e errors in line fluxes 


(0.01-0.05)H/3 


30.5 


42.5 


(0.05-0. 15)H/3 


22.5. 


34.5 


(0.15-0.30)H;g 


15.5 


27.5 


(0.30-2.0)H,9 


10.5 


16 


(2.0-5.0)H;g 


8.5 


10 


(5.0-10.0)H/3 


6.5 


8.5 


> lOH/3 


5 


8.0 




Electron Densities 


(cm ^) and Temperatures (K) 


Ne[S II] 


600±90 


600±140 


Ne[Cl III] 




450±270 


Te [O III] 


9100±2960 


11550±4950 


T,. [N 11] 


S500±1950 


9550±3340 


Te [S 11] 


9150±1250 


8400±3010 


" In units of 10~~ 


^'^crg cm^"^ s~^ 





The electron temperature in the knots is measured in this 
paper for the first time. 

Summarising, the densities and temperatures of the two 
knots are the same within the errors. 



6 NGC 7662 

This bright nebula (PN G106. 5-17.6) consist of a central 
cavity surround ed by three co ncentric shells, the inner Rim, 
the outer shell ([Guerrero. Jax on & Chu 2004), and a halo 
(jCorradi et alj [20041 ') . In addition, on smaller scales it has 
a number of LIS. Most of them appear, in projection, dis- 
tributed at the outer edge of the outer shell. All these struc- 
tu res are easily ident ified in Figure El an d are named as 
in lBalick etHf (119981 ) and IPerinotto et all (i2004,). This fig- 
ure also identify the structures th at we analyze in this pa- 
per. iBalick. Preston fc Ickd (|l987l ) studied the kinematics of 
the different nebular components, from which the acronyms 
FLI ERs and SLOWERs w ere later introduced. 

IPerinotto et al.l (|2004h studied this PN with unprece- 
dented spatial resolution using the STIS spectrograph on 
the HST. However, their spectra were taken through two 
slits (called POSl and POS2, see their Figure 1), that in 
terms of LIS, are restricted to only one side of the nebula. 
We have used much longer slits that cover all the structures 
along both sides of the nebula, in each direction at P.A. 175° 
and 248° . This allows us to compare the densities and tem- 
peratures of any given symmetrical pair of LIS. The slit at 
P.A. 248° contains a pair of low-ioniz ation filaments (th e 
serpentine-like filaments described in iBalick et al.l 119981 ). 
The slit at P.A. 175° does not contain any LIS, and is used to 
compare the physical parameters along the line joining the 
FLIERs with those in a direction free of micro structures. 
In this way, the p resent work can be see n as complementary 
to the analysis of IPerinotto et al] (|2004| ). 

In Table 5 we report the derived values of ch/3 that are 
based on the average weighted Ch/3 given by the Ha, H7, 
and H5 to H/3 ratios. These values are all between 0.13 and 
0.21 for the 5 regions under analysis in each position angle. 
These ch/j agree very well with those previously published : 
0.18 - 0.22 (iTvlenda et al]|l992l ): 0.1 (iHvung fc Alledll997l ): 
0.2 (|Perinotto et al. "2004'): and 0.18 lZhang et al.l (12004 ^ 

This table gives the observed line fluxes, densities and 
temperatures of the 10 structures we measured along P.A. 
175° and 248°. These parameters are also presented in Fig- 
ure (2] The average Te [O ill] derived from the emission in- 
tegrated along both P.A. are the same within errors (Te[0 
lll] = 14 300±700 K and 13 400±680 K, respectively). In fact 
these values also agree with the Te[N n] = 12 750±1 530 K 
derived along the P. A. =248°. The Te[0 ill] wc obtained here 
are more in agreement with the results by Barker (1986) - 
whose [O III] temperatures in different regions of the nebula 
(see his Tab le 3) v aries from 11 200 K to 13 800 K - and 
I Zhang et all (|2004|) fl3 300 K ) , tha n with the somewhat 
lower values of IPerinotto et al] (|2004l ) . 

As for the [S 11] and [CI lll] densities of the integrated 
emission we found 2 300±440 cm"^ and 2 550±490 cm"^ 
for P.A. =175° and 2 850±450 cm"^ and 2 850±460 cm"^ 
for P. A. =24 8°. NefS 11] results very closely reproduce those 
found by a IWang et atl (|2004 Ne[S ll]=2884 cm'^'). The 
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Figure 5. NGC 7662: HST archive images. Top: [O III]. Bottom: 
[N ll]. In both images the box corresponds to 34x34 arcsec^. 
The size of the extraction windows in the spectra are as follows. 
P. A. =175° (inner rim + outer FLIER): -14'.'7 to -lO'.'S, SE outer 
shell; -8'.'4 to -2'.'8, SE inner rim; -1'.'75 to +1'.'75, Cavity + Star; 
+2'.'8 to +S'.'4, NW inner rim; +12'.'6 to +14'.'7, NW outer shell. 
P. A. =248° (inner rim + outer shell): -14'.'7 to -11'.'2, SW FLIER; 
-9'.'8 to -2^8, SW inner rim; -2^1 to +2^1, Cavity + Star; +3'.' 5 
to +9^8, NE inner rim; +10"5 to +15'.'4, NE FLIERS. 



same authors give a value of 1 862 cm~^ for the [CI ill] 
density, but errors are not reported in their work. 

Focusing on the densities and temperatures of the inner 
Rim and outer shell with respect to those of the FLIERs 
(pair of filaments at P.A.=248°), it is clear from Figure [21 
that: i) temperatures Te[0 ill] (where errors are smaller) are 
slightly higher at the inner Rim than in the outer knots 
(P.A.=175°), but this is just at the 2-a level of the adopted 
errors; ii) densities of the Rim agree very well one to another 
in both position angles; iii) the Rim is roughly a factor of 2.8 
denser than the outer shell (P.A.=175°); iv) the Rim (inner 
shell) along P.A.=248°, is a factor of ~ 1.7 denser than its 
outer region and t he pair of FLIER s. 

bummansmg, IPerinotto et~aLl l|2004l ) previously ob- 



served the LIS (FLIERs and SLOWERS) of this PN at 
different directions than in the present work. Densities for 
the FLIERs in both works are in good agree ment: Ne=2 
200 c m~^ for the two FLIERs stud ied bv iPerinotto et al] 
(|2004 ). and 2 100 cm~^ in our work. IPerinotto et al.l (|2004h 
did find significantly different densities for FLIERs and 
SLOWERs, being the latter approximately 1.3 times denser 
than the former. 

Our analysis shows that FLIERS of NGC 7662 are less 
dense than the inner regions of the nebula, but they are at 
the same time much denser than the PN outer shell in which 
they are located. 



7 DISCUSSION 

7.1 Temperature of the knots 

We have estimated individual temperatures Te[0 ill] and 
Te [N ll] for the different regions of the four nebulae studied 
in this paper. As it is well known, temperature errors are 
always important and they range in our case typically from 
around 1000 K for Te[0 ill] in the brighter zones to more 
than 4000 K for Te[N ll] in the fainter regions. Therefore, a 
sensible comparison of temperatures for the different com- 
ponents of each PN is precluded, and only rather general 
conclusions can be extracted: 1) temperatures are those of 
normal photoionized PNe, varying from 9 500 to 14 500 K 
typically; 2) Te[0 III] agrees within the errors with Te[N ll] 
for every PN component measured. 3) No relevant differ- 
ences between Te[0 ill] at the LIS and at the higher excita- 
tion PN structures are found. 



7.2 Knots density contrasts 

The empirical density contrasts between the knots and the 
inner nebular regions, as described in this paper, are such 
that knots are typically 1 to 3 times less dense than the rims 
or cores of these PNe (He 1-1, IC 2149 and NGC 7662): in 
no one of the LIS measured here the electronic density is 
definitely larger than at the main PN structures. This un- 
expected result is neverthel ess also found in t he ful l sample 
of PN with LIS. Table 4 of IConcalves et all (|200ll ) hsts 26 
PNe with pairs of knots similar to those in the present work. 
We have searched in the literature for density measures both 
at the inner rims (or shells) of each PNe and at their pairs 
of knots, finding eigh10 PNe (in addition to the four neb- 
ulae here studied) with adequate data. In all but one PN 
(NGC 6826) the densities of the knots are equal or lower 
than those of the inner rims or shells, by factors between 1 
and 10. 

This result is not consistent with the theoretical ex- 
pectations. The formation of PN major structures and 
pairs of jets/knots has been modeled in the past using 
a broa d variety of processes and scenarios (see, for i n- 
stance. IConcalves et al] (120011 ) and lBahck fc FrankI l|2002l )'). 



NGC 4634 ICuerrero. Jaxon fc Chul |2004|); Hb 4 
Haiian et al]|l997l); I C 4593 llCorradi et al.lll9 97'): NGC 7009 
Goncalves et al.'l2003!); K 3-35 iMiranda et al.''200 0); NG C 6826 

(Balick et al. 1994) ; N GC 24 40 (.Cucsta fc Phillipa,2000li ; K 1-2 

fexter, PoUacco fc Bellll2003ll . 
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Table 5. NGC 7662 Observed line fluxes, flux errors, cng, Ne and Tg. P.A.=175° (Rim and 
Outer shell) and P.A.=248° (Rim + FLIERs) 



P.A.=175° Line ID 


SE outer Shell 


SE inner Rim 


Cavity + Star 


NW inner Rim 


NW outer Shell 


NEB 


H5 4101.8 


30.81 


27.31 


28.38 


29.48 


24.25 


28.99 


H7 4340.5 


23.79 


23.99 


22.00 


23.99 


24.60 


23.99 


[Olll] 4363.2 


17.17 


15.66 


19.78 


18.46 


19.12 


17.82 


H/3 4861.3 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


[Olll] 4958.9 


504.4 


368.0 


343.2 


338.7 


481.7 


372.1 


[Olll] 5006.86 


1533. 


1105. 


1018. 


1018. 


1434. 


1114. 


[CLill] 5517.7 


0.761 


0.317 


1.756 


0.345 


0.398 


0.424 


[CLill] 6537.9 


0.662 


0.344 


1.669 


0.378 


0.475 


0.423 


[Nil] 5754.6 


0.169 


0.087 


0.469 


0.099 


0.152 


- 


Ha 6562.8 


331.6 


324.9 


334.9 


326.7 


334.9 


339.9 


[Nil] 6583.4 


7.496 


2.981 


2.502 


2.640 


5.191 


3.116 


[Sll] 6716.6 


1.164 


0.366 


0.242 


0.249 


0.916 


0.371 


[Sll] 6730.8 


1.428 


0.693 


0.373 


0.396 


1.107 


0.636 




3.164 


36.90 


14.15 


30.94 


2.138 


104.4 


CH/3 


0.18±0.017 


0.16±0.012 


0.20±0.016 


0.17±0.015 


0.18±0.017 


0.21±0.02 








Percentage errors in line fluxes 






(().()l-0.(J5)H/3 


35 


15.5 


12.5 


16.5 


30 


8.6 


(0.05-0. 15)H/3 


16 


10.5 


10 


9.5 


11 


4.6 


(0.15-0.30)H/3 


10 


8.6 


9 


8.5 


9 


4 


(0.30-2.0)H/3 


7 


6.6 


7.6 


6 


7 


3.6 


(2.0-5.0)H/3 


5.5 


4.6 


6 


5.5 


6.6 


3.6 


(5.0-10.0)H/3 


5 


4 


4.6 


5 


6 


3.6 


> lOH/3 


4.5 


4 


4 


4.5 


6 


3.6 






Electronic Densities (cm 


and Temperat 


ures (K) 




Ne[S 11] 


130t)±735 


3650±1()60 


2950±730 


3400±1()50 


1200±560 


2300±440 


Ne[Cl III] 


1250±7()() 


3450 ±1000 


2100±370 


3550±1100 


1100±550 


2550±490 


Te[0 III] 


1240()±135() 


1345()±122() 


15600±1600 


15000±1390 


13200±1300 


14300±700 


Te[N 11] 


12650±4830 


13900±3350 




16850±4250 


14450it5300 




P.A.=248° Line ID 


SW outer 


SW inner 


Cavity + Star 


NE inner 


NE outer 


NEB 




FLIER 


Rim 




Rim 


FLIER 




[Sll] 4068.6 


1.420 


0.885 


0.814 


0.920 


1.263 


1.232 


[Sll] 4076.4 


0.630 


0.469 


0.966 


0.430 


0.317 


0.269 


[Sll] 4072.0 


1.923 


0.964 


2.127 


1.131 


1.381 


1.804 


H6 4101.8 


24.39 


24.98 


25.08 


23.69 


23.93 


23.99 


H7 4340.5 


43.96 


46.22 


45.07 


44.17 


43.99 


46.29 


[Olll] 4363.2 


21.83 


17.27 


16.43 


16.53 


18.89 


16.98 


H/3 4861.3 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


[Olll] 4958.9 


552.7 


387.6 


339.8 


398.1 


517.8 


397.1 


[Olll] 5006.86 


1706. 


1184. 


1040. 


1214. 


1570. 


1213. 


[CLill] 5517.7 


0.637 


0.395 


0.408 


0.389 


0.606 


0.410 


[CLill] 6637.9 


0.686 


0.431 


0.390 


0.437 


0.570 


0.422 


[Nil] 5754.6 


0.685 


0.088 




0.089 


0.552 


0.166 


Ho 6562.8 


343.8 


331.9 


321.6 


331.3 


333.8 


321.1 


[Nil] 6583.4 


31.14 


3.661 


2.734 


4.289 


22.88 


6.321 


[Sll] 6716.6 


2.630 


0.464 


0.311 


0.431 


2.141 


0.611 


[Sll] 6730.8 


3.843 


0.761 


0.423 


0.679 


3.064 


0.932 




3.580 


33.69 


22.23 


41.22 


8.732 


117.9 




0.21±0.026 


0.17±0.011 


0.14±0.011 


0.17±0.011 


0.19±0.015 


0.14±0.008 








Percentage errors in line fluxes 






(0.01-0.05)H;3 


38.5 


10.5 


18 


9.5 


26.6 


8.6 


(0.05-0. 15)H;3 


19.5 


8.6 


9.5 


7 


14 


6.6 


(0.15-0.30)H/3 


13 


6.6 


7.0 


5.5 


8 


4.6 


(0.30-2.0)H/3 


8.5 


4.6 


6.0 


4.5 


6 


3.6 


(2.0-5.0)H/3 


6.5 


4.0 


4.6 


4.0 


4.6 


3.0 


(6.0-10.0)H/3 


5.0 


3.6 


4.0 


3.5 


4 


3.0 


> lOH/3 


4.5 


3.6 


3.6 


3.5 


4 


3.0 






Electronic Densities (cm 


and Temperat 


ures (K) 




Ne[S II] 


2400±1300 


3950±650 


1850±550 


3250±550 


2250±800 


2850±450 


Ne[Cl III] 


1800±1040 


3500±580 


2150±630 


3800±1370 


2000±790 


2850±460 


Te[0 III] 


14750±1980 


13650±970 


14100±1070 


13300±820 


12700±1100 


13400±680 


Te[N 11] 


11350±4540 


12600±1920 




11700±3920 


12900±3720 


12750±1530 



In units of 10 



Note that [Nil] 6548A was not measured for any of the structures. Thus, Tg [N ll] 



were calculated assuming the theoretical relation between [Nil] 6583A and [Nil] 6548A. 
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Figure 6. Diagnostic diagrams for spatially resolved PNe. Lef t panels are an ad aptation of the lRaga et al. I l l200^ results, whereas right 
panels superpose to these results the data for our sample. The lRaea et al.l ||2008| ) simulations for high- velocity shocked knots that travel 
through a photoionized region of a PN involve four families of models, as indicated in the labels of the upper left panel. Coded in 
their names, the numbers represent the knot's velocity, from 40 to 150 kms~^whereas letters from 'a' to 'C indicate different stellar 
luminosities and temperatures which are kept constant in each model while the knot's distance to the ionizing source decreases. Therefore, 
model 'alOO' has the lowest, and 'ClOO' the highest, photoionization ra te, while both assume a velocity of 100 kms~^for the knots. A 
compilation of empirical literature data performed bv lRaga et al. ] lj2p08i) is also sketched in the four panels as regions filled with diagonal 
lines. In the two right panels, different regions of the PNe of our sample are added: empty symbols correspond to the pairs of knots, 
filled symbols to the nebular inner regions (Rim, Core, Cavity + Star), and half-filled circles correspond to the outer shells. 



Recent models can be grouped into three main families. 

i) Magnetohydrodynamic disk wind models, based on the 
scenajio of accretion disks formed by binary interaction s 
l|Frank fc Blackmanll200i . iBlackman. Frank fc Weichll200lf ) . 

ii) Interacting AGB and post-AGB wind models, in which 
the post-AGB wind is driven exclusively by the magnetic 
press ure of a single star (Ga rci'a-Segura, Lopez & Franco 
2005. iGarci'a-Di'az et aUbOOSl ). iii) Models that explore the 
effect of axisymmetrical hght (low density) jets that, for a 
short period of time, expand into a spherical AGB wind 
l|Akashi fc Sokeill2008l ). Some of the above models are par- 
ticularly tailored to reproduce the main morphological and 
kinematic properti es of a few well observed P Ne: He 3-401, 
M2-9 and He 2- 90 (iGarci'a-Segura et aLlbOOSl ): Mz 3 and the 
pre-PN Ml-92 (|Akashi fc SokeJl2008h . Nevertheless, and de- 
spite the different assumptions of the different models, all of 
them are able to account for the formation of PNe with jets. 

Model predictions state that, after the eventual shut off 
of the jets, dense knots would form embedded in the outer 
nebular components. Moreover, both the knots present at 
the tips of the existing jets and those left over when jets dis- 
appear are predicted to be denser than the inner PN compo- 



nents by large factors ranging from 10 up to 1000 (typically 
10) depending on the model assumptions. 

Therefore, there is a clear disagreement between ob- 
served and modeled pairs of knots: the former are found 
to be typically 10 times less dense, while the latter are ex- 
pected to be typically 10 times denser, than the inner rims 
and shells of their parent PN. 

A straightforward possibility to reconcile the empirical 
and theoretical density contrasts is keeping in mind that 
empirical densities, as those measured here, correspond to 
only the ionized fraction of the PN gas. Models, on the other 
hand, usually predict the formation of knots accounting for 
its morphology, kinematics and total gas content, but do 
not detail the ionization status of each component nor their 
evolution with time. The strong discrepancy between mod- 
els and observations found here (a typical factor of 100) can 
therefore be explained if 90% or more of the low-ionization 
knots' material is neutral. It is worth mentioning here the 
case of one prominent cometary globule of the Helix neb- 
ula -which contains the best known low-io nization knots 
of PN e. The molecular matt er (H2 -|- CO: iHuggins et al] 
l2002h and the dust content (|Meaburn et al.lll992l 'l of this 
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knot amount to ~2xlO^^M0, whereas lO'Dell fc HandronI 
1 19961 ) estimated that the typical mass of the photoionized 
gas of a globule in Helix is ~1O~^M0. 

In another wel l-studied system of knot s, NGC 7009, pre- 
vious modeling by iGoncalves et al] (|2006l ) succeeded in re- 
producing the observed ionization structure and chemical 
abundances under the assumption that the total gas density 
of the knots is equal to the one em pirically derived from the 
[S ii] lines (jCongalves et al.ll2003h . therefore, not account- 
ing for any neutral matter. This apparently contradicts the 
present results suggesting that the knots might be mostly 
neutral. However, the effe ct of the density assumed for the 
knots in the modeling by IGoncalves et all l|2006h is partly 
balanced by its geometry, i.e., by its shape and size. As the 
shape was fixed based on the knots' appearance in the HST 
images, the size was allowed to vary until it fits the observa- 
tions (i.e. t he optical line ra t ios an d electron density). For 
that reason IConyalves et al.l ()2006h were able to reproduce 
the observations of NGC 7009 without invoking need for 
neutral matter. It would be ver y interesting to a t temp t 3D 
photoionization modeling as in IGoncalves et all l|2006t) ex- 
ploring high density regimes for the knots as suggested in 
this paper. 

7.3 Main excitation of the pairs of knots 

iRaga et al.l l|2008l ) have recently simulated the evolution of 
high-density (10^ cm"^), high-velocity (100 to 150 kms"^) 
knots that travel away from a photoionizing source crossing 
a uniform lower density (10^ cm~'^) environment. They ex- 
plored a range of initial conditions for the shocked knots (in 
terms of photoinization rate), and qualitatively reproduced 
the emission-line ratios of several FLIERs, rims and shells 
observed in pre-PNe and PNe. At the end of the model evo- 
lution, 400 yr, the number densities for the knots are 10^ to 
10'' cm~^, broadly in agreement with measured Ne of knots 
in NGC 7009 and IC 4634. These simulations were indeed 
intended to reproduce the line ratios of the FLIERs of these 
two PNe. No predictions are given for the properties of the 
other PN components (rims, shells and haloes). In Figure 6 
we sh ow two of the diagnostic diagrams from iRaga et al.l 
(|2008t ) comparing empirical and theoretical line ratios for 
FLIERs, as well as showing a compilation of empirical data 
for rims and shells of some PNe. 

Fig. 6 shows, first, that most of the inner regions (rims 
and shells) of the PNe in our sample occupy the same zone 
as the rims and shells of other PNe in the literature. Sec- 
ond, the same is true for the zone occupied by the pairs of 
knots in He 1-1 and NGC 7662 (open boxes and circles) with 
respec t to the zones of the FLIERs sample in iRaga et al.l 
l|2008t ). Third, knots in the other two nebulae, IC 2149 and 
KjPn 8 (open triangles and stars), are d isplaced from the 
empirical zone of FLIERs in iRaga et all (2008). Their lo- 
cations are consistent with the results of simulations with 
the highest initial velocity of the knots (150 kms~^, model 
A/B/C/150), and indicate that a mixing of shock and photo 
ionization is playing a role in the measured emission-line ra- 
tios. In fact, t he highly superson ic velocity (radial velocity 
of 226kms~^ iLopez et all 119971 ) measured in KjPn 8 for 
the SE knot at P.A.=98° could explain this agreement. Ex- 
pansion ve locities of IC 2149's knots are more moderate (~ 
47 kms~^ lvTzauez et al.ll2002l ). 



iRaga et al.l l|2008l ) show that models with low pho- 
toionization rates yield line ratios typical of shock-excited 
knots, whereas in models with higher photoionization rate, 
photoionized emission dominates the spectrum of the high- 
velocity, high-density kno ts. This strengthen the idea that it 
is evolution (PN age ; see lGoncalvesll2004 iRaga et al."2008|: 
Viironen et al.ll2009l ) and photoionization rate (|popita 1997; 
Raga et al. 20081 ) . rather than the highly supersonic veloci- 
ties of the knots, the key parameters to separate the shock- 
excited from the photoionized structures in planetary nebu- 
lae. 



8 CONCLUSIONS 

We have derived the physical parameters of low-ionization 
pairs of knots, and of higher ionization PN main components 
for a small sample of nebulae (He 1-1, IC 2149, KjPn 8 and 
NGC 7662). We have shown that electron temperatures at 
the knots are comparable to the temperatures of the main 
nebular components, rims and shells, whereas densities in 
the latter structures are significantly higher than in the pairs 
of knots, by up to a factor of 2. Typical knots' densities in 
the studied PNe are 500 to 2 000 cm"^. 

We have argued that optical line ratios are not appropri- 
ate to constrain the total gas density, as they only account 
for the ionized matter of the knots, whereas none of the 
available models properly explore the knots ionization frac- 
tion, but the total (neutral plus ionized) density structure. 
A direct comparison of our results, plus compiled literature 
data, with the theoretical density contrasts expected from 
model predictions, yields a clear discrepancy that amount, 
typically, up to a factor of 100. To reconcile models and 
observations we suggest that an important fraction of the 
gas at the knots -90%- should be neutral. The results pre- 
sented here can help in constraining more realistic models 
for jets and knots in PNe, by providing empirical physical 
parameters of the ionized gas. 

We have analyzed the location of th e knots of our four 
PNe in diagnostic diagrams studied by iRaga et al.l (j2008l ) 
concluding that the observed line emission ratios are com- 
patible with shocked knots traveling in the photoionized 
outer regions of their host PNe. 
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